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Abstract. The basic physical assumptions and results of Landau's hydro- 
dynamic model of particle production are reviewed. It is argued that these 
results have sufficient descriptive and predictive power in strong-interaction 
phenomenology, including recent RHIC data, to warrant a closer examination 
of the physical assumptions. 
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1. Introduction 

The success of boost- invariant hydrodynamics at RHIC (e.g. [Q]) to describe the 
systematics of elliptic flow is considered to be strong evidence that the strongly 
interacting system thermalizes early, coupling the dynamical evolution to initial- 
state geometry of the system. These calculations can, in principle, provide a means 
to extract the form of the equation of state of the thermalized system. This may 
provide insight into whether or not a quark-gluon plasma was in fact created in 
heavy ion collisions at RHIC. 

It should not be forgotten, however, that hydrodynamics has a long history 
in the study of strongly-interacting systems. Coupled with the intense interest in 
statistical and thermal model calculations in the early 1950's, spearheaded by Fermi 
and Landau [El El 0, this led to a large body of strong interaction phenomenology, 
manifestly non boost-invariant, which was refined throughout the years [El E| 
IHl El ^1 O- It is the goal of this talk to review the central physical assumptions 
and predictions of the Landau hydrodynamical model (and its refinements), and 
to show its relevance to a wide range of results in a variety of strongly-interacting 
systems. 

In the Landau-Fermi physical picture, the main physical assumptions were: 



• The collision of two Lorcntz-contracted hadrons or nuclei leads to full thermal- 
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ization in a volume of size Vrup/y^. This justifies the use of thermodynamics 
and estabhshes the system size and energy dependence. 

• A massless blackbody EOS is assumed p = e/3. This allows the complete 
calculation of physical quantities. 

• All chemical potentials (meson and baryon) are zero, which dramatically sim- 
plifies the mathematics. 

The main results derived from these assumptions are: 

• A universal formula for the produced entropy, determined mainly by the initial 
Lorentz contraction. 

• Gaussian rapidity distributions, at least for particles produced several units 
away from the projectile rapidities 

• Thermal particle occupations determined by T ~ m^r . 
2. Universal Entropy 

Simply using the first law of thermodynamics and the blackbody EOS, Landau and 
Fermi both arrived at the same scaling formula for the multiplicity produced in a 
collision of two strongly interacting objects, Nch = aS — Ks^l'^ [El El 12 ■ The 
agreement between this formula and a wide range of systems is shown in Fig. ^ 
with a constant of proportionality of K=2.2, which has been found to work with 
systems as diverse as p+p, Au+Au, e+e^ and + p [0 El Ej. 

It was always controversial to apply Landau's reasoning to "small" systems, 
such as e~ annihilations. Beyond the usual arguments about falling far short 
of the thermodynamic limit, it has been argued that perturbative QCD (pQCD) 
already provides an excellent description of not just the multiplicity of produced 
hadrons as a function of energy, but the details of jet shapes (modulo hadronization 
corrections), rendering a "statistical" description superfluous (at best!) Despite 
this, it is interesting to observe that pQCD calculations (e.g. Ref. [E|) and the 
Fermi-Landau formula agree within 10% over a wide range of energies, essentially 
up to where the LEP data ends, as shown in Fig. ^ (adapted from Ref. [Ej)- This 
similarity is striking when one considers that pQCD implies an infinite mean-free- 
path for parton rescatterings, while the hydro description implies a negligible one 
for the (presumably partonic) degrees of freedom which thcrmalize. 

The extension of the Fermi-Landau approach from a single p+p collision to nu- 
clei is surprisingly simple. Since the Lorentz contraction is not changed, the angular 
distributions are in principle similar to the smaller system [0}. Then, provided that 
the interactions between the subvolumes of the system do not themselves open up 
new degrees of freedom, the entropy of the system as a whole will simply scale pro- 
portionally to the number of participating nucleons (Npart)- This "A'part-scaling" 
has been observed for the total multiplicity (but not necessarily in any particular 
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Fig. 1. Charged particle multiplicities for A+A, p+p (with leading particle effect 
removed) and e+e^. Theoretical curves are pQCD (dotted line) [E], baryon-free 
Landau (solid line) [ E| , and Landau including the baryochemical potential effect 
(dot-dashed) 



region of phase space) in all collisions involving nuclei, from p + A to Au + Au 
collisions [[HI [HI- 

3. Thermal Phenomenology and Hadrochemistry 

In the Landau scenario, freezeout is not expected to occur immediately, as Fermi 
assumed, but rather when the temperature reaches the limit of the pion Compton 
wavelength T = rriT^. This was based on a suggestion by Pomeranchuk [ 115) to avoid 
Fermi's prediction that nucleons would outnumber pions by virtue of their larger 
statistical weight. This assumption leads to predictions for the relative population 
of various particle states similar to those made in the Hagedorn approach [ I16L I17 |. 

A+A collisions clearly deviate from the Fermi-Landau formula at low energies. 
An obvious suspect is the phenomenon of baryon stopping, which is absent in p+p 
collisions but is substantial in A+A [El- If one puts back the —^bNe term into the 
first law of thermodynamics, we immediately see how the presence of a conserved 
quantity associated with a substantial mass (i.e. the proton mass) will naturally 
suppress the total entropy: S = {E + pV — iibNb)/T. Using an existing thermal 
model code, Cleymans and Stankiewicz [ I19| calculated the entropy density as a 
function of ^/s. It rises to limiting value where fiB ^ and T ^ Tq, the Hagedorn 
temperature. If we then assume that the total multiplicity scales linearly with the 
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total entropy, then the suppression relative to the asymptotic value should describe 
the energy dependence of the A+A multiplicities relative to the Landau formula. 
This allows a direct comparison with the multiplicity data as shown in Fig. ^ 
There we find a reasonable qualitative description of the difference between A+A 
and e^e~ collisions. 



4. Gaussian Rapidity Distributions 

Up to this point, the discussions have only involved the total entropy, integrated 
over the full phase space. Landau was apple to perform approximate calculations of 
the angular distributions by introducing hydrodynamic evolution using the standard 
equations of relativistic hydrodynamics d^T^^" = closed by the blackbody EOS 
[13 El- These equations generically imply that the initial state entropy, which is 
produced in the process of thermalization, is distributed in rapidity space with a a 
Gaussian form, the width determined by the initial Lorentz contraction. Including 
the total multiplicity formula to set the overall normalization, the full expression [ 
El is: 

dN Ks^/^ , y2 

where L = ay = (1/2) ln(s/77ip) = ln(7). Already, one can see a connection between 
the shape of the distribution and the multiplicity, since from this definition s^/^ cx 

While the experimental status of Landau Gaussians vs. Feynman-Bjorken 
plateaus [1201 1^ in dN/dy was ambiguous throughout the 1970's [El El) despite 
strong evidence for the applicability of Landau's formulas, the situation in A+A 
collisions was clarified rather quickly by the two large-acceptance RHIC experi- 
ments, PHOBOS and BRAHMS. PHOBOS quickly estabhshed that boost invari- 
ance was violated over a large rapidity range by inclusive measurements of dN / drj 
over jyyj < 5.4 [,22,. BRAHMS consolidated these observations by finding that 
the rapidity distributions of pions in |?;| < 3 is Gaussian with a width parameter 
only 10% different from the Landau prediction [ES|- This led to a re-evaluation of 
all of the existing heavy ion data, where it was found that charged pion rapidity 
distributions all fell close to the Landau trend [ I24L I25| . 



5. Connections or Coincidences 



While the Landau results seem describe, and even predict, several non-trivial results 
at RHIC, the formulas also encode features that are usually attributed to QCD, or 
its various approximations. 

For example, measurements of dN/drj in p+p and A+A collisions, boosted 
into the rest frame of one of the projectiles, showed the phenomenon of "limiting 
fragmentation" [ 22 . This is energy- independence of the particle yields at a fixed 
rapidity distance from the projectile rapidity [ I26| . One thing that has not typically 
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Fig. 2. BRAHMS data on dN/dy 
for charged pions, compared with Lan- 
dau's prediction. Inset is a comparison 
between the measured and predicted 
Gaussian widths (from Ref. [I23p 
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Fig. 3. Resuhs from the Landau for- 
mula, Equation in the rest frame 
of one of the projectiles. Approximate 
limiting fragmentation is observed. 



been appreciated in the context of the Landau approach is that limiting fragmenta- 
tion seems to arise naturally from Equation : Transforming this distribution into 
the rest frame of one of the projectiles, y' — y — ybeam, one finds that approximately: 

Which only varies weakly with L for y' ^ 0. This seems too strong to be merely 
a coincidence, but it is not clear why limiting fragmentation, which implies xp 
scaling, arises naturally from the Landau formulas, in which xp does not seem to 
be a preferred variable [ |H] . 

Another interesting, and unexpected, connection can be made between the Lan- 
dau expressions and the results from the CGC-based model of Kharzeev and Levin 
[ I27| . Once one fixes the peak dN/dy of both models, and then varies the energy, 
it appears that the width of dN/dy varies in a similar way as a function of beam 
energy, as shown in Fig. ^ It is interesting that the exponent extracted from 
HERA data, A — 0.25 — 0.3 is surprisingly similar to the power seen in the Landau 
multiplicity formula. However, it is not clear why this similarity occurs, or how 
robust it is. 

A final unexpected coincidence is seen in the transverse direction near y=0. 
Carruthers and Duong- van noticed that the pr distribution of tt^'s in p+p collisions 
was well described out to pT = 10 GeV by a Gaussian distribution in transverse 
rapidity yr — \ l'^( ™T-pr ^ ^ ^ ^'^^ [EHl- While no derivation was given 

for this phenomenological description, which holds over 10 orders of magnitude. 
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Fig. 4. Comparison of the Landau 
result, Equation^ compared with CGC 
results from Ref. [EZI- The results are 
normalized at y=0 at 200 GeV. 
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Fig. 5. Top: STAR data on i^d^A^/dp^ 
and PHENIX data on EdPa/dp^ fit to a 
Gaussian in yT- Bottom: Data divided 
by the fit. 



an argument was made on a similar basis as for the Gaussian in the longitudinal 
direction. To see if this function continues to work well at RHIC energies, fits have 
been made to PHENIX 7r° data [1^ and STAR inclusive charged data from 
p + p collisions, which are shown in Fig. |S1 Reasonable agreement is found with the 
STAR data with L = 0.56, despite the combination of various particle species, and 
excellent agreement is found with the PHENIX data with L = 0.54, up to pT = 11 
GeV. 

6. Conclusions 

In conclusion, it is argued that the arguments made by Landau and collaborators 
in the mid-1950's appear to have a surprising relevance for understanding RHIC 
phenomena. The results on the total multiplicity and the shape of the rapidity 
distributions hold in a robust way for a variety of systems, including p-\-p, e+e" and 
Au+Au, perhaps providing a natural way to understand the apparent universality 
seen in the total number of charged particles. Several interesting connections are 
found between the Landau results and limiting fragmentation, CGC calculations, 
and even particle production at very high pT- Understanding these connections may 
provide deeper insight into the strong interaction and the dynamical properties of 
strongly interacting systems. 
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